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Abstract. We compute the transmission of an electron through an impurity in 
polyacene. For simplicity the disorder is confined to a single unit cell. When the 
impurity preserves the inversion symmetry around the central axis, the scattering 
problem can be reduced to that of two independent chains with an alternating sequence 
of two types of atoms. An analytical expression for the transmission coefficient is 
derived. On-site and off-diagonal defects are considered and shown to display very 
different electron scattering properties. 
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1. Introduction 

The study of polyacene started long ago in a pioneering paper by Kivelson and 
Chapman. pQ The study performed by these two authors was motivated by the 
investigations of Su, Schreiffer, and Heeger on polyacetylene. [2J [3] TheypQ showed that 
the band structure associated to the 7r-orbitals of polyacene is made of four bands, with 
a valence and a conduction bands touching at the limit of the one- dimensional Brillouin 
zone. Since polyacene can be considered the extreme case of a narrow graphene [H [5] 
ribbon it is instructive to compare the band structure of graphene and polyacene. In 
both cases one has an electron per carbon atom and in both cases the valence and 
conduction band touch at the corner of the Brillouin zone. In graphene, however, the 
bands are linear at the K and K' points of the Brillouin zone, whereas in polyacene 
the dispersion relation around the momentum k = n/a (a the length of the unit cell) 
is quadratic. As a consequence, graphene has a vanishing density of states at the 
Fermi energy whereas polyacene displays a square root singularity. It was therefore 
speculated [TJ that the ground state of polyacene could support long range order such as 
ferromagnetic and superconducting phases. It was further showed that phonon modes 
in polyacene can lead to a Peierls distortion opening an energy gap at k = ir/a. [HE] 
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Using a Green's function method, Rosa and Melo[7] studied the density of states 
of several distorted configurations of polyacene, including the case where a local defect 
was present. They found that the optical response of polyacene should be very different 
from that of polyacetylene. The inclusion of many-body effects on the description of 
the ground state of polyacene was done by Garcia-Bach et al. [8] using a valence-bond 
treatment. They found that the distorted ground state is degenerate. Using a projector 
quantum Monte Carlo method, Srinivasan and Ramasesha p] found, within the Hubbard 
model, that electron-electron correlations tend to enhance the Peierls instability. With 
the development of the density matrix renormalization group, the exact study of the 
ground state of quasi-one- dimensional electronic systems became available. Raghu et 
al. [10] studied the ground state of polyacene using the Pariser-Parr-Pople Hamiltonian. 
As in previous studies they found that strong electron-electron interactions can enhance 
the dimerization. Using a configuration interaction technique, Sony and Shukla [TT] 
studied the optical absorption and the excited states of polyacenes also in the framework 
of the Pariser-Parr-Pople Hamiltonian. 

As interesting as the ground state nature of carbon polymers are their transport 
properties. Of particular interest to us is the effect of disorder on the electronic tunneling 
through a portion of a disordered chain. Guinea and Verges, [12] using a Green's 
function method, studied the local density of states and the localization length of a 
one-dimensional chain coupled to small pieces of laterally linked polymer. They showed 
that at the band center there is an exact vanishing of the transmission coefficient due 
to a local antiresonance. Sautet and Joachim [13] studied, within a one- dimensional 
tight-binding model, the effect of a single impurity which would change both the on-site 
energy and the hopping to the next-neighbor atoms. Mizes and Conwell [14] consider 
the effect of a single impurity in the square polymer, showing that a change on the 
on-site energy has a more pronounced effect in reducing the transmission in the one- 
dimensional chain than in the square polymer. At the end of the paper [2] these authors 
speculate that for polyacene there should be four active scattering modes instead of two 
as in the square polymer. As we show in this paper, this is not the case because the 
band structure of polyacene is markedly different from that of the square polymer. Gu 
et al. [15] generalized the study of Ref. [14] by including different types of obstacles as 
scattering centers. Yu el al. [16] studied the electronic transmission through a conjugate- 
oligomer using both mode matching and Green's function methods. They found that the 
transmission through the conjugate-oligomer possesses several transmission resonances. 
Farchioni et al. [17] studied the transport properties of emeraldine slats using a Green's 
function method, by looking at the effect on the electronic transmission of impurity 
dimers. The band structure of very narrow carbon nanoribbons was studied by Ezawa. 
[j~8] The effect on the conductance of the parity of the number of carbon rows across a 
zig-zag ribbon has been studied by Akhmerov et aL[T9] 

In this paper we study the effect of disorder on the electronic transmission in 
polyacene. This system can be considered the most extreme limit of narrow carbon 
nanoribbons. [T8] and it corresponds to an odd parity situation studied in Ref. [19J. 
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We shall assume that the disorder is limited to a single unit cell. Although this 
assumption can be relaxed, it allows us however to obtain a full analytical expression 
for the transmission coefficient. Following Ref. [13] we shall consider both on-site and 
hopping disorder. The system is therefore characterized by two semi-infinite perfect 
leads made of polyacene and a scattering region. The Schrodinger equation has to be 
solved for the three regions and a mode matching technique, first used by Ando [20] 
to study the conductance of a square lattice in a magnetic field, will be used here. 
Although, as discussed above, all the theoretical studies point out that polyacene has 
a gap at the Fermi energy due to phonons and interactions, here we will consider non- 
distorted polyacene in the independent particle approximation. It should be relatively 
simple to generalize the calculations below to include the effect of a gap and the effect 
of electron-electron interactions (at least at the mean field level), but the only expected 
modification would be a small region near zero energy where the conductance would be 
zero due to the presence of the gap. In addition, the disorder could add states into the 
gap of polyacene. 

2. Model for disordered polyacene 

The model for disordered polyacene can be written in general terms as 

oo 4 

H = J2 [$ijt( n , i) \n, i) (n, j | 

n=— oo i,j=l 

+ ti ntIa ,(n,i,j)\n,i){n,j\ 

+ t intcv (n,i,j)\n,i)(n + l,j\] , (1) 

where e(n,i) represents the on-site energy, ti QtTa (n,i, j) the hopping between the atoms 
within the unit cell n, and ^^(71,1, j) the hopping between atoms in neighboring unit 
cells. The sum over n runs over the unit cells and the sums over % and j run over the 
atoms in the unit cell. In our model we shall assume that the on-site energies e(0, 1) = eo 
and e(0,4) = eo are different from those of the rest of the atoms in the lattice. The 
on-site energies for atoms to the left of these two are all equal to ei and, for those to the 
right, to €r. This difference could be due to a potential bias applied to the system. The 
hopping parameters are defined in Figs. [D and El Only the hopping energies starting 
at the atoms 1 and 4 of the zero unit cell are modified. They are named ti and t^. 
We further assume that the hopping energies between atoms 2 and 3 within a given 
unit cell (t±) can be different from that between atoms 1 and 2, or 3 and 4 (t). The 
eigenfunctions of this problem have the form 

I^) = EE^,j)Kj)- (2) 

n 3 

As it is clear from Fig [2] the problem is divided into three regions: (i) that to the left 
of the zero unit cell, with n = — oo,...,— 1, representing a semi-infinite ordered lead; 

(ii) the disordered unit cell at n = 0, which will induce the scattering of the electrons; 

(iii) the right ordered lead, for n — 1, . . . , oo. We have therefore to solve the problem in 
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(a) 



(b) 



Figure 1. (color online) Representation of polyacene in (a) and its topologically 
equivalent lattice in (b). Each unit cell has four atoms. The on-site energy is cl.r- 
The hoping along the chain is t and in the perpendicular direction is t±. 




Figure 2. (color online) Characterization of the n — unit cell, where the disorder 
is located. Note that the on-site energies for n < (el) are different from those for 
n > (er). Also the on-site energy of atoms 1 and 4 in the zero unit cell is different 
from the rest and equals eg. 



these three regions. The method we will use was developed by Ando[20j in the context 
of quantum point contacts in a magnetic field and generalized by Khomyakov et al. [21] 
to an arbitrary three-dimensional structure. 

The strategy of solution is therefore the following, firstly the problem is solved 
for the pure leads obtaining the eigenmodes. Then the scattering problem is solved 
considering an incoming wave described by one of the eigenmodes at a time. This 
allows for the determination of the transmission matrix elements t^ u , and from them 
the calculation of the conductance follows as [22, 23J 




where the summation is over propagating channels only. 
3. Solution of the scattering problem 

As explained above, the problem is separated into the solution in the perfect leads and 
the solution in the scattering region. In what follows we shall explain how to obtain 
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those solutions. 



3.1. Scattering channels 



We study here the scattering channels, i.e. stationary solutions in the asympotic leads 
assuming these are perfect and infinite. Let us start with the left lead. The problem for 
the right one is solved along the same lines but with replaced by sr. The tight-binding 
Hamiltonian for m = — oo, . . . , —2 has the form 

- -Bc m _i + (IE - H L )c m - B^c m+ x = , (4) 

where c m is a 4-vector containing the wave-function coefficients of the unit cell m. The 
matrices B and are given by 
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It is important to note that the B matrix is singular, not having an inverse. This will 
be important in what follows. As explained by Ando, [20] the problem for the perfect 
lead may be solved assuming a Bloch relation between the vectors c m 



(7) 



For a general problem the vector c m has dimension M and the solution of (j4]) can 
be obtained by transforming it into an ordinary eigenvalue equation of dimension 
2M.[20j We look for scattering states formed by an incoming wave approaching the 
scattering region from the left and its resulting outgoing waves. Whether these waves are 
propagating or evanescent states depends on the nature of the eigenvalues. Propagating 
states always have |A| = 1, whereas evanescent ones have |A| ^ 1. Another possibility is 
to transform Eq. (jlj) into a quadratic eigenvalue equation by repeatedly use of Eq. ([7]), 
resulting in 

-Bc m + (IE-H L )\c m -B^\ 2 c m = 0. (8) 

As before, the determination of the eigenvalues of Eq. (JSj) will lead, for a problem of 
dimension M, to a polynomial of order 2M whose roots are the sought eigenvalues. 

Although we can attack the solution of the problem ([8]) using the matrices given by 
Eqs. ([5]) and (jSJ), it is however convenient to perform a unitary transformation of Eq. 
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). We define the unitary matrix U 

I 1 1 
J_ 1 1 

V2 



(9) 



10 0-1 

V o i -l o / 

and introduce the transformation c m = Uc m , and M = UMU^ 1 . Upon this 
transformation, the eigenproblem (jBJ) is reduced to two block-diagonal quadratic 
eigenvalue problems, since the resulting transforms of B and are factorized into two 
2x2 matrices which we denote generally as b and h^. Correspondingly, c m becomes the 
direct sum of two 2-vectors which we generically refer to as u m and which we assume 
normalized. The resulting eigenvalue problem reads 

-bu m + \{IE-h L )u m -\ 2 tfu m = 0, (10) 



with 



and 




:ir, 



JJ- (12) 

The decoupled two-dimensional problems describe propagation through the even and 
odd modes with respect to the central axis of the polyacene. The eigenvalue problem ( JTOl) 
has the same form as that of a linear chain, with hopping energy t between neighboring 
atoms and with on-site energies alternating between and ± t±. As a conclusion, 
the scattering problem in polyacene can be mapped into that of two independent one- 
dimensional chains, contrary to the expectations of Mizes and Conwell. jH] 

From the general discussion above one would expect that the eigenvalue problem 
(II Op would lead to a quartic polynomial in A. In fact because the matrix B has no 
inverse (as happens for b) the polynomial is only cubic in A. One of the solutions is 
the trivial one A = 0. This solution has to be disregarded since it would produce a null 
wave function everywhere. The other two solutions are 

A = [{E - e L f + t ± (E - e L ) - 2t 2 ]/2t 2 
1 



± 2 V K E " ^ + ~ ^) ~ 2t2 ] 2 A 4 - 4 ■ (13) 
When the square root becomes imaginary, Eq. (fT3|) gives the momentum of a 
propagating Bloch wave. If we now consider the case t± —>■ —t± [see Eq. ( fT2l) ] two 
other solutions are obtained. 

The fact that we only have two solutions for each sign of t±, and not four, means 
that there are always two of the four expected modes that do not contribute to the 
transport, not even as evanescent waves. This result could have been anticipated if we 
had considered the energy bands of perfect polyacene. pQ These are given by 

t± , 1 



E - e L = ±-± ± -yjl? ± + 8t 2 [l + cos(fca)] , (14) 
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where a is the length of the unit cell and k G [— n/a, it /a}. If we now try to solve for k 
in Eq. ffT4l) one finds that for a given energy E only two bands give a solution, being it 
real or complex. 

The velocity of the electrons in the modes is given by [21] 
2a 

For the present problem the velocity has a simple form given by 



[Au t 6 t 



u 



(15) 



(16) 



where ua and ub are the components of u m , orbital A resulting from the (symmetric 
or antisymmetric) linear combination of sites 2 and 3 within the cell and orbital B 
stemming from the similar combination of sites 1 and 4 (see Fig. [I]). These amplitudes 
are given by 

t\l + \\ 



(E-e L ) 2 + t 2 \l + \\ 



and 



u B 



(l + \*)(E-e L ) 
|l + A| v /(E-e i ) 2 + i 2 |l + A| 2 
Using Eqs. (TTTj) and fTT8l) the velocity ffT6l) reads 
2at t(E - e L ) 



-3A. 



h (E-e L ) 2 + t 2 \l + \\ 



(17) 



;is) 



(19) 



Equation (Tl9|) allows to identify the right and left moving modes for a given energy E. 
The right lead is solved in the same way with replaced by €r. 



3.2. The scattering region 

We now want to describe the scattering region. The Schrodinger equation for the unit 
cell m = — 1 has the same form as before, except that it couples to Cq. For the unit cell 
m = the Schrodinger equation is written as 

- Be , + (IE - flooH - B ] rCi = (20) 

with the matrices Hqq and Br given by 
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For the unit cell m = 1 the Schrodinger equation has the same form as Eq. (J3J) except 
that B is replaced by B R , is replaced by e R and it couples to c . As before we can 
perform a unitary transformation of the Schrodinger equation leading to an effective 
2x2 Hamiltonian. The Hamiltonian for unit cell m = — 1 has the same form as Eq. 
(|T0l) . For the unit cell m = we obtain 

- bu ! + (IE - h 00 )u - b RUl = , (23) 

and for m — 1 

- b R u + {IE - h R )ui - tfu 2 = . (24) 

The matrix b R has the same form as b with t replaced by t R . The matrix h R is obtained 
from Hl replacing by e R . The matrix h 00 is given by 

,X)- (25) 

Since the full problem factorizes into two block-diagonal problems and, furthermore, 
because in the leads, for each 2x2 block, only one mode is active, the scattering 
takes place without mode mixing. This is a considerable simplification over the general 
approach of Refs. [HEt]. 

If we define as the amplitude at cell m of the Bloch wave propagating to the 
right (+) or left (— ) in perfect polyacene, the following relation holds: 



u 



m+l 



A ± ^ m , (26) 

where A 1 * 1 = 1/A T is given in Eq. ( |T3l) for the asympotic left lead. Using Eq. ( 1261) we 
can write 

u 2 = (A" - A + ) u±! + A+it_i . (27) 
The boundary conditions require the specification of ul^ which will represent an 
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Figure 3. (color on line) Representation of the transmission T{E) for different values 
of en- The parameters are t = t± = tL = tR = —2.7 eV, eo = e L = 0. The solid line 
corresponds to perfect polyacene. 

incoming wave function in one of the modes of the left lead. On the right lead we have 
only a scattered wave propagating to the right, we thus write 

u 2 = X + Ul . (28) 
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Figure 4. (color on line) Representation of the transmission T(E) for different values 
of eo- The parameters are t = t± = ti, = in = —2.7 eV, £l = £r = 0. The solid line 
corresponds to perfect polyacene. 



Using Eqs. fl27l) and fl28|) . the determination of the wave function on the unit cells 
m — — 1, 0, 1 reduces to the resolution of the following system of linear equations 

[IE -h L - A + b]u_i - tfu = {\- - X + )ut 1 , (29) 

- bu ! + (IE - h 00 )u - b ] RUl = , (30) 

- b R uo + (IE - h R - 6 f A + )wi = . (31) 

Once the system is solved the vector Ui is determined. We can then write U\ as 

u x = w+r Q , (32) 

where a is an explicit index to label the (odd or even) mode. The physical transmission 
matrix element is given by[23j 

ta = \ Ta i (33) 

V v L 

where vl/r represents the velocity of the considered mode in the left /right lead, given 
by Eq. (fl9|) . The total conductance is obtained by summing the contributions from the 
symmetric and the antisymmetric mode, 

2p 2 2p 2 2 

G = \T ( E)^ K \>, (34) 

a=l 

assuming that both modes are propagating waves. Solving explicitly the linear system 
of equations defined above, the expression for r a is obtained, 

r a = qtH L t R u\X~ — , (35) 

where 

q = \~ - A + , (36) 
E — e R 

V = — , 37 

uZ 
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and the m-independent parts of the A and B components of u+ in the left 

and right leads, respectively. The denominator D reads 

D = - t\\t 2 + {t 2 - E{E + U) + e L {2E + t ± - e L ))\ L ] 
x [t 2 - E(E + t ± ) + e R (2E + t ± - e R ) + t 2 \+] 
+ {E + t ± - e L ) [-t 2 + {-2t 2 + E(E + *,_) - e L {2E + t ± - e L ))X L ] 
x [t 2 R {-E + e R ) + (E — e )(-t 2 + E(E + t ± ) - e R {2E + t ± - e R ) - t 2 \$$) 

We are now in position to study the transmission coefficient T(E). In Figure [3] we 
study the transmission of the electrons in a situation that mimics that of a potential step 
in ordinary quantum mechanics problems. The potential barrier is created by having 
the sites at the right of the unit cell n = at a different energy from those at the 
left. For definiteness we consider = 0. When = e R , T(E) has a step structure, 
represented by the solid line in Fig. [3J due to the existence of two possible conducting 
channels. These two conducting channels are those associated with the two effective 
one-dimensional chains. For energies where the two channels of the two chains are 
propagating states, one obtains T(E) = 2; when only one channel is active one obtains 
T(E) = 1. Having now ^ e R induces some back scattering of the electrons at the 
interface, reducing the value of T(E). When the difference between ei and e R increases 
we see the appearance of zones of zero transmission, while this is always nonzero for 
£l = e R . 



t R =-3.7eV 

t R =-4.7eV 

trf=-1.7eV 

t R =-0.7eV 



■"''/ r/ Pv 



-5 5 

energy £ (eV) 

Figure 5. (color on line) Representation of the transmission T{E) for different values 
of the hopping and ti,. The parameters are t = t± = —2.7 eV, £l = £r = to = 0, 
tft = t^. The solid line corresponds to perfect polyacene. 

In Figure H] we study the effect of changing the value of eo relatively to that of 
and e R . Overall there is an effect of diminishing T(E). If the difference between eo and 
€l and e R is not very large the T(E) curve follows closely that for the non-disordered 
case. There is however an exception at zero energy, where T(E) tends to zero, leading to 
a totally reflecting barrier. This can be due to a local antiresonance as in Refs. [T2l [24] . 

On the contrary, if the disorder is induced by changing the hopping, the behavior 
of T(E) shows almost perfect transmission at E = 0, as can be seen in Fig. [5j This 
effect is particularly clear when ti and t R , become very small. From Eq. (1351) we see 
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that T(E) will vanish as the fourth power of ti (ti = nevertheless close to E = 
the system has a strong enhancement of its transmission. Away from this special point 
the curves for T(E) are similar to those of Fig. H]for on-site disorder. 

In the cases of Figs. 0] and [5] the conductance of the system would exhibit 
conductance oscillations as long as the scattering strength of the defect is not too strong. 

4. Conclusions 

In this paper we have studied analytically the electronic transmission of polyacene due 
to on-site and off-diagonal disorder. We have shown that the system is equivalent to 
two decoupled linear chains of atoms, with alternating on-site energies and a constant 
hopping parameter. Scattering occurs in the system without mode mixing if inversion 
symmetry around the central axis is preserved. We find that on-site and hopping 
disorder have markedly distinct scattering properties close to zero energy. Whereas 
for on-site disorder the transmission decreases strongly close to zero energy, for hopping 
disorder the transmission is enhanced. If the system opens up a gap at the Fermi energy, 
we expect these same characteristics to occur near the gap edge. 
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